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Radiation of a charged particle moving parallel to a inhomogeneous surface is considered. Within 
a single formalism periodic and random gratings are examined. For the periodically inhomogeneous 
surface we derive new expressions for the dispersion relation and the spectral-angular intensity. In 
particular, for a given observation direction two wavelengths are emitted instead of one wavelength 
of the standard Smith-Purcell effect. For a rough surface we show that the main contribution to 
the radiation intensity is given by surface polaritons induced on the interface between two media. 
These polaritons are multiply scattered on the roughness of surface and convert into real photons. 
The spectral-angular intensity is calculated and its dependence on different parameters is revealed. 

PACS numbers: 41.60.-m,42.25.Fx,42.79Dj,41.75.Ht 



Introduction. Smith-Purcell radiation(SP) [1| is origi- 
nated when a charged particle travels parallel to a plane 
with diffraction grating. Recent renewed interest in this 
problem is caused by different applications. Among these 
applications are length determination for short electron 
bunches [3, creation of monochromatic light source in 
the far infrared region [3-13 1 ^^'^■- Various theoretical 
models were proposed for describing the S P: [9Hl3| , for a 
brief review of recent theoretical works see [Ij, [1^1 ■ Most 
of these models deal with the periodical grating in the 
strong scattering regime (see below). However in many 
situations the interface over which the charge travels is 
rough . As an example one can mention chamber walls 
in storage rings. Even the best treated surfaces contain 
roughness. Radiation appearing when a charged particle 
moves near a rough surface could be useful for beam diag- 
nostics [16] . The influence of the surface roughness on the 
transition radiation (originating when particle crosses the 
interface between two media) was discussed in [13, [3 ■ 
Roughness-induced radiation for a charged particle slid- 
ing over a surface was experimentally observed in [19| . 
In the present paper we study radiation emitted due to 
electromagnetic field scattering on inhomogeneities of di- 
electric constant . We will see below that in the weak 
scattering regime it is possible to develop a rigorous the- 
ory describing both periodical and random grating within 
a single formalism. 




FIG. 1: Geometry of the problem. A charged particle moves 
parallel to the Oa; axis. Observation point is far away from 
the system. 



Maxwell equation for the electric field reads 

V^E(f,io) — gTa,dd[vE(f,uj) + -^e(r,uj)E('r,Lu) = j(r,uj) 

c^ 

(1) 

where j is the current density related to the charge 



j{r,.)^-'-^5[z-d)8[y)e-^l^ 



(2) 



General Relations The geometry of the problem is 
shown in Fig.l. 



A charged particle moves uniformly in the vacuum at 
the distance d from the plane z — Q separating vacuum 
and isotropic medium. We are interested in the radiation 
field far away from the charge and the interface. The 



Here v is the velocity of the particle moving on Oa; 
direction and e(r, w) is the inhomogeneous dielectric per- 
mittivity of the system which for a rough surface can be 
chosen in the form 

£(r, w) = 9[z — h[x,y)\ + e[uj)Q[h{x,y) — z\. (3) 

where 9(z) is Heaviside's unit step function, h{x,y) is 
the amplitude of surface roughness. As it follows from 



Eq.dS]), the space z > h{x,y) is vacuum while the space 
z < h{x, y) is occupied by a medium with isotropic di- 
electric constant £{uj). Assuming h{x,y) small and ex- 
panding Eq.dS]), one gets [20| 



e{r,uj) =£n{z,u) + £r{r,uj) 



where 



and 



£q(z,uj) 



1, z>0 

£{uj), z <Q 



£r[r,uj) = [£{uj) - l]S{z)h{x,y). 



(4) 
(5) 
(6) 



Thus the total e is presented as a sum of a regular part £q 
and an irregular part Sr- To separate the radiation field 
we decompose electric field E = Eq + Er analogous to 
Eq. dH) . Here Eq and Er are the background and radiation 
fields, respectively. They obey the following equations 



V EQ{r,u!) — graddivi?o(r, w) - 
+^r£o{z,uj)Eo{r,Uj) = j(r,w) 



(7) 



V'^Er(f,io) - gra,ddiv Er(f,uj) + -^eo(z,uj)Er(f,Lu) + 

& 

+^y£r{r,uj)Er{r,^) = :^£r{r,uj)Eo{f,uj) (8) 

c c 

Note that although the term £rErin Eq.® is small one 
should keep it because it causes multiple scattering of 
electromagnetic field. We will see that multiple scatter- 
ing effects are very important in radiation from rough 
surface. Multiple scattering effects in SP radiation for a 
cluster of dielectric particles were discussed in [2l|. At 
large distances from the system the electromagnetic field 
can be treated as a plane wave in which electric and mag- 
netic fields equal to each other. Therefore the intensity 
of radiation at the frequencies [uj,uj + dcu] and at solid 
angles [n,n + dH] can be determined as follows 



dl{uj,n) = -\Er{R)?R^dQ.duj, 



(9) 



where n is unit vector on the direction of observation 
point R, fi is the corresponding solid angle, see Fig.l 
and also [1J|. As usual at large distances \Ey.{R)\'^ be- 
haves as \/R? therefore intensity does not depend on 
R. The expression Eq.® should be averaged over the 
realizations of random roughness h{x,y). For this rea- 
son it is convenient to introduce the Green's functions of 
Eqs.dMD 

£o{z,u})^-dxf_c ~ -T, — ^— + 

+Sx^.y'] G%{r,/,u) ^ SxJ{r- r^) (10) 

£o{z,uj) — dxf^ - ^^ ^^ -f- oa^' 



drxdvu 



Sxu^^- 



+£r{r,uJ)—^x^, 



G/iJ.(f,r',w)=,5A,<5(r~r').(ll) 



In Eas. (llOlll|) a summation over the repeated indexes 
is supposed. Solutions of inhomogeneous Eas. (|7l8p can 
be expressed through the Green's functions Eqs. (|10|lip . 
Using Eqs. (|7l8l) and (jlOllip one can represent the 
averaged radiation intensity tensor < Iij(R) >=< 
E„{R)E*^{R) > in the form 



< I.jiR) >^-r drdr' < G,^(i?,r)e,(f) 

G:^ir^,R)£rin > Eo,if)E*,{r') (12) 



where the background electric field Eof^{r) is expressed 
through the bare Green's function 



Eo^if)^ dnGlxir,n)jxiri) 



(13) 



Here < . . . > means averaging over the surface random 
profile h{x,y). Note that in the original Smith-Purcell 
experiment [1], as well as in subsequent works on SP a 
periodical grating in one direction is used. In this case 
h(x,y) = h[x) is some periodical function of one coor- 
dinate. In the present paper within a single approach 
we consider both periodical and random gratings. In the 
random case we suppose that ft. is a gaussian stochastic 
process characterized by two parameters 



< h{p) >= 
<h{pi)h{p2)>=5^W{\pi-p2\) 



(14) 



where p is the two dimensional vector in the xy plane, 
5"^ =< ft^(/o) > is the average deviation of surface from 
the plane z = 0. Correlation function W is character- 
ized by a correlation length a at which it is essentially 
decreased. 



The Maxwell equations for electric fields Eq. (|7l8p 
and Green's functions Eq. pOllip should be amended by 
boundary conditions. As usual, it is required that tan- 
gential components of electric field be continuous across 
the plane z = 0. The exact field, of course, will satisfy the 
boundary conditions across the surface z = h{x,y) rather 
than the plane. However this approximation seems rea- 
sonable for small roughness A 3> (5 and is widely used in 
the literature. The Green's function 6*^1^(^,7^,0;), when 
considered a function of z for fixed z' satisfies the same 
boundary condition as the pth Cartesian component of 
electric field. 

Green's Functions. The equation for bare Green's func- 
tion Eq. ljlOp with correct boundary conditions for arbi- 
trary £{ui) was solved in ^20]. To obtain radiation inten- 
sity in vacuum we will need Green's functions in the half 
space z > 0. In order to simplify the problem we will 
consider the case when isotropic medium is a metal with 
very large negative dielectric constant |e(w)| ^ 1. Using 
expressions for Green's functions from [20| we find the 



following basic components 



iri^ p(i,Ap^1^ 



°"<*'> -="(*•»' -1^^.), 



G',M'-.ii) = -GUm--} = 



iqz 



Wx £[uj)qe 
k^ ki — e{(jj)q 



(15) 



where G^Ap\z,z') is the two-dimensional Fourier trans- 
form of G^j{r,r') and z > 0. In the coordinate represen- 
tation 



G?-(f,r') = 



dp 

f2^ 



,w(P-ff)Gl{^z,z') (16) 



Here p and p are two-dimensional vectors with Carte- 
sian components px,Py,0 and x,y,0. Also k = uj/c, fci 
and q are determined as follows: 




fci 



p^, k^ > p^ 
- k^, k^ Kp^ 

-(e(^)fc2_p2)i/2 



(17) 
(18) 



In Eg. p^ a branch cut for the square root along the 
negative real axis is assumed [20|. Other components 
of Green's function are small over the parameter 1/|£|. 
To determine radiation intensity we will need asymp- 
totics of Green's functions at large distances. Substi- 
tuting Eg. lfTS)) into Eg.dTSl). one finds 



G^,(i?,p,0) 



1 



2ttV2R 
1 



(k{R-npp) + ^) 

y/rT^ cos I k{R — Upp) + — 



gIMp,o) = -gUp,o,R)~ 

If / ^ ^ TT 

y^sin ( fc(i? — fipp} + — 



2ttV2R 
rix 



\k(R - Upp) - - 



■ sm 



{k(R 



- rioP - 



2ttV2R 



flpTlxTlz 



sin [k{R - ripp) - ^) - 



^-^ sin ( k(R — Hpf 
/n7 V 



Radiation Intensity. Spectral-angular radiation inten- 
sity Eg . p^ can be represented as a sum of two contribu- 
tions, I{R,uj) — I°{R,uj)+I^{R,uj), where /q and /^i are 
sing le scattering and diffusive contributions, respectively 
[231 ■ First consider the single scattering contribution to 
the radiation intensity. Substituting the Green's func- 
tions in Eq. ([T2l) by the bare ones, we obtain 



(19) 



where n is the unit vector on the direction of the obser- 
vation point R = fiR, nx,z and Up = 



corresponding components. When obtaining Eq.(|19p we 
use asymptotics of Bessel functions for large argument 
[2^ . Eqs.([T9l) are correct provided that kR:^\, Rp:^ p 
and we use approximate equation |i? — r| ss R— fir. 



I^^iR) = {e- ifS'k^ J dpdf^GliR, p, 0) 
Gf^iff, 0, R)W{\p- ff\)Eo4p, 0)E*,iff, 0) (20) 

where (ij) = (xz). The background electric field in the 
limit |e(a;)| > 1 can be found from Eqs.(I13|), ^ and JH]) 



EoMO) = - 



4ggifeoa; 



dkc) 



V 7vV+^ 



^^^WZ±Z) (21) 



7 



where fcg = oj/v, 7 = (1 — v'^/c^)^^^^ is the Lorentz 
factor of the particle and A'l is the first order Mac- 
donald function. As it follows from Eg. ipTj) the back- 
ground electric field and correspondingly radiation in- 
tensity is exponentially small when cod/vj ^ 1, see also 
|25| . One can expect essential intensity provided that 
ujd/vj ^ 1. Far away from the system at the ob- 
servation point one can use asymptotic expressions for 
Green's functions Eg. p^ . Substituting Egs. ([T5|) and 
(pri) into Eg. ipO)) . for the spectral-angular radiation in- 
tensity I{uj,n) = cR^Iii{R)/2, one obtains 



/"(w,rj) 



e^ gLx{l^nl){l + nl){l+nl) 
c/32 



IfSiYUpd 



(22) 



where L^ is the system size in the x direction, g — 
{e - l)25V2fc4 and P = v/c. When obtaining Eq.(|22|) 
we neglect strongly oscillating terms in the limit kR ^ 1 
and suppose that W{p — ff) = (j'^5{p — f/). Beside that 
we substitute the Macdonald function by it's asymptotics 
for small argument assuming that kod/j <§; 1. In the op- 
posite limit as was mentioned above radiation intensity 
is negligible. The components of unit vector n are de- 
termined through the polar 9 and aziniuthal angles 
of observation direction: n^ — cos 0, Tip — sm9,nx — 
sin 6 sin 0. We consider radiation into the half-space 
z > (vacuum) which means 9 < tt/2. Note that the 
couphng constant g = fc'*(e — l)^S^a^ in Eq. f^ is a di- 
mensionless parameter. From the condition Rp 3> p, 
one obtains a restriction on angles sm9 ^ L/R, where 
L is a characteristic size of the system. To avoid mis- 
understanding note that 1//?^ dependence of radiation 
intensity Eq. (j22p is correct in an intermediate regime for 
not very low velocities uid/vy ^ 1. When /3 — ?> 0,as was 
mentioned above, radiation disappears. 

Note that the background field Eg can originate ra- 
diation without any roughness provided that Chcrenkov 
condition v^e > (? is fulfilled. Cherenkov radiation is 
possible for dielectric surfaces with positive large e. For 
metallic surfaces in the optical region we are interested 



in the present paper dielectric constant is negative and 
Cherenkov radiation is absent. 

Periodical case. Analogously one can consider the case 
when surface grating is a periodical function. For sim- 
plicity we will assume that h{p) = Sshi2TTx/b, where 
b is the period of grating. Substituting W{\p — p*|) by 
sm2'Kx/bsm2TTx' /b in Ea. ipO]) and using Eg. ((TO)) , after 
integration , for spectral-angular radiation intensity, one 
has 



Isp{n,uj) = 
5{knx + fco 



e^ g,{l + nl){l-nl){l+nl)L., 
- -r) +(>{ko - kn^ r^ 





Fikny) (23) 



where gi — k (e — l) S and F is determined as follows 

2 



Fikny) 



dkn 



1 Jo 



dy 



Kii 



feovz+z^ 



Akn-i^ 



ViFTlP 



(24) 



When obtaining Ea. (j23|) we keep only the terms propor- 
tional to Lx- In the most interesting case Uy ^ and 
dfco/7 *^ 1: substituting Ki by its asymptotic expression, 
one finds from Eq.(IM|),F(0) ^ n^/A. As follows from 
Eg.ip^. because of the 6 functions, for a given observa- 
tion direction, only two discrete wavelengths are emitted 



A± = b{- ±11^) 



(25) 



This is a generalization of well known Smith-Purcell dis- 
persion relation (ll| to the weak scattering (see below) 
case. For an arbitrary periodical grating one can ex- 
pand the surface profile h{x) into Fourier series and 
for each term one can obtain analogous dispersion re- 
lation with b substituted by b/m, where m is the diffrac- 
tion order. Note that the dispersion relation Eg. ([25]) 
and the spectral-angular radiation intensity Eg. ip^ dif- 
fer from reported earlier. The reason of those differences 
are following. First we are considering weak scattering 
regime instead of strong one considered in above men- 
tioned papers. Our theory is applicable provided that 
{e{uj) - l)'^5^/\'^ < 1 although |£(a;)| » 1. Probably 
this regime was realized in the experiment on SP radi- 
ation in the optical region for shallow gratings [2J]. As 
mentioned in [24] traditional formula of SP radiation is 
failed to explain the results of experiment in the shal- 
low grating case. Second reason is the boundary condi- 
tions. As follows from Eqs. (l7l8p Maxwell equations for 
Eq and Er contain the same disruptive function £0(2;, 1^). 
Therefore both of them should satisfy the same bound- 
ary conditions at z = 0. In our consideration this goal 
is achieved automatically because the Green's functions 
Egs. (|15ll9p satisfy the correct boundary conditions |20| . 
In contrary, in traditional consideration [111 ], only the 
total field Ef) + Er and not they separately, satisfy the 
boundary conditions at 2; = 0. Probably this difference 
leads to different dispersion relation Eq. 



P2 



PmnlSPl P2) 



FIG. 2: Dashed line is the correlation function of roughness 
g5{pi — P2) and the solid line is the averaged over the random- 
ness two-dimensional Green's function of surface polariton 



Diffusive Contribution, Surface Polaritons. Using 
Eg. ([12)) one finds diffusive contribution to the radiation 
intensity in the form 

ij^ (R) ^ g dpidp2dpGi,n {R, Pi , 0)Ghz {p2, p) 

Pmnhsipi - P2)G:M P2)G:^(0, pi , i?) |i?0. (P, 0)\^ (26) 

where Gij{p2,pi) = G'ij(/>2, O,/?!, 0),and where diffusive 
propagator P is determined by the sum of ladder dia- 
grams see Fig. 2. and 26]. 

All integrations over z coordinates make them equal 
to because of 5{z) in fluctuation part of dielectric con- 
stant Eq.®. Averaged two-dimensional surface polari- 
ton Green's function [23] satisfies the Dyson equation 



G^u{p)^GlM+9Glr.{p) 



7^~^G^n{pi)Gniy{p) 

(27) 
Remind that Gp,v{p) = G'^i,(p]0, 0),see Eg. (|T5)) . Bare 
Green's functions are determined by Eg. (ITS)) . In further 
we will interested in the behavior of Green's function 
close to the pole. These values give the main contribu- 
tion in the limit g — )■ 0. As it follows from Eq. lfTS)) two- 
dimensional Green's functions of surface polariton has a 
pole at p^ = fc^e/(e -|- 1), see [23]. The corresponding ve- 
locity of a surface polariton is equal to c-\/(£ + l)/£ < c. 
Remind that we consider the case when £ ^ — 1. When 
electron velocity becomes equal to this velocity a super- 
radiant emission is possible provided that the grating is 
periodical [y, Q . Close to the pole and for large negative 
]£(a;)] ^ 1 the Green's functions of surface polariton can 
be represented in the form 



gIAp) 



gUp) 



-k 



1 



v/=i7H^'-p' 



la 



Wx 



^-£i(cj) k"^ -p^ - ic 



(28) 



where a = k e2/e{, e = £i+i£2 and £2 ^ l£il. InEq.(| 
a describes the damping of the surface polariton on the 
flat surface due to the inelastic processes in the medium, 
i.e. £2(0;). It follows from Eq.dMj) that / dppG°^ (p) = 0. 
Therefore only G^., gives contribution to the integral in 



Ea. (l27l) . Solving Dyson equation Ea. (l27|) one can repre- 
sent the averaged Green's functions in the form 



gUp) 



-k 



1 



iPx 1 



(29) 



where A = / -J^ImGl^{p) == gfc^/4£i(i:j). Real part 
of the integral leads to renormalization of the parameters 
and does not play any role. Integral is calculated in the 
limit £2 — >■ 0. A describes the damping of the surface 
polariton by its roughness-induced conversion into radia- 
tive modes [27|. It is convenient also to introduce the 
polariton mean free path on the rough surface I — k/A. 
Note that the neglected terms in diagram expansion are 
small on parameter \/l <^ 1 [26| . 

Using Eas. ((29|) one sees that the main contribution to 
the diffusive radiation intensity Eg. d^Bl) give the term 
proportional to Pzzzz which contains a diffusive pole at 
small momentums. Summing the ladder diagrams in 
Fig. 2 one finds a Bethe-Salpeter equation for diffusive 
propagator P{K) = Pzzzz{K) 



P{K)^g^f{K)+g!{K)P{K) 



where 



I{K) 



dp 



■G{p)G*{\p~K\) 



(30) 



(31) 



Here G{p) = Gzz{p)- Using Eq. l^^ and calculating the 
integral in Eq. (|31]) in the limit g — >■ 0, one finds P{K) at 
small Kl <t:l 



P{K) 



2g 



(32) 



Substituting Eqs. ([T9| into Ea. (l26l) . for the diffusive con- 
tribution to the radiation intensity, one has 

1° (w, n) = — ^ :^^^ ^P{K -^ 0) 



GATT^Tln 



dp\Eoz{p,0)\' 



(33) 



As follows from Eas. ((32|) and (p3)) radiation intensity 
diverges. This divergence is caused by the infinite size of 
the system, see also [2j, [2^ . If one takes into account the 
finite size of the system the minimal momentum will be 
of order ^ 1/L. As was mentioned above the radiation 
intensity is exponentially small provided that d 3> 7/fco- 
In the opposite limit d ^ 7/fco substituting Ki by it's 
asymptotic expression and integrating Ea. (|55)) . we finally 
obtain 

^ (.'^'"''-7^ ^ZZ IJT- '.34J 



c/32 



8TTn„ 



dp 



In this consideration weak I <C ?m, where U 



k/a 



£i/ke2 is the inelastic mean free path of surface po- 
lariton, absorption can be taken into account as fol- 
lows [29]. When L > (Ihnf^, L in Eq.dM]) should 
be substituted by {llinY''^- Comparing single scatter- 
ing Eq. ((22| and diffusive Ea. (p4l) contributions, one has 
/^//" -- L^/P > 1. Therefore the diffusion of surface 
polaritons is the main mechanism of radiation. Let us 
make some numerical estimates for the optical region. 
For Ag at \^ 4500A° £1 7.5 and £2 - 0.24. Tak- 
ing 5 - mA° and a - 1000A° [2J7,i, one has g ~ 0.68, 
I ^ 7.7A and Un ^ 47.94A. Thus the conditions for dif- 
fusive mechanism X <ti I <ti lin, L are fulfilled. Evidently, 
depending on grating parameters S, a emission in other 
wavelength regions is possible too. 

We have considered multiple scattering effects in radi- 
ation for uncorrelated roughness. However they are very 
important for the periodical as well as correlated grating 
cases too. These cases are more complicate and will be 
discussed elsewhere later. Our result Eq. ([23| for SP ra- 
diation intensity with only single scattering contribution 
is correct in the cases when multiple scattering contri- 
bution is negligible. Such a situation can occur for the 
metals with relatively large absorption when the condi- 
tion of multiple scattering of polaritons Im 3> I is not 
fulfilled. 
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